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Abstract. The NLO QCD calculation for the photoproduction of the isolated 
photon with a large pT at the HERA ep collider is presented. The single resolved photon 
contribution and the QCD corrections of order as to the Born term are consistently 
included. The sizeable NNLO contributions, the box and the double resolved photon 
subprocesses, are taken into account in addition. The importance of the isolation cut, 
as well as the influence of other experimental cuts on the px and rj^ distributions 
are discussed in detail. Results are compared with experimental data and with the 
different NLO calculations. 
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1. Introduction 

The production in the ep colhsion of the prompt photon with large transverse momentum 
Pt is considered. Such reaction is dominated by events with almost real photons 
mediating the ep interaction, ~ 0, so in practice we deal with the photoproduction of 
the prompt photon (called also Deep Inelastic Compton (DIG) scattering). The photon 
emitted by the electron may interact with the proton partons directly or as a resolved 
one. Analogously, the observed final photon may arise directly from hard partonic 
subprocesses or from fragmentation processes, where a quark or a gluon decays into 7. 

The importance of the DIG process in the ep collision for testing the Parton Model 
and then the Quantum Ghromo dynamics was studied previously by many authors [|l|- 



Measurements were performed at the HERA ep collider by the ZEUS group [|T0 



1^ , and |jT3|, also the HI GoUaboration has presented preliminary results fl^. In 



these experiments events with isolated photons were included in the analysis, i.e. 
with a restriction imposed on the hadronic energy detected close to the photon. The 
corresponding cross sections for the photoproduction of an isolated photon and of an 
isolated photon plus jet were calculated in QGD in next-to-leading order (NLO) [pTSf - 
2^ . There exists analogous calculation for the large-Q^ case (DIS events) p3[ . 

In this talk the results of the NLO QGD calculation for the DIG process with an 



isolated photon at the HERA ep collider are presented [jT8|, ^ . We consider the parton 
distributions in the photon and parton fragmentation into the photon as quantities of 
order aem- We emphasize the importance of the box diagram jg —>■ jg, being the 
higher order process, in description of the data. Our approach differs from the NLO 
approach |2l], ^ by set of subprocesses included in the analysis. The comparison 

of our predictions (KZ) |]l9l with the results obtained by L.E. Gordon (LG) [jl^ and M. 
Fontannaz et al. (FGH) and with data measured by ZEUS group is presented. 

2. The NLO calculation for jp — > jX — general discussion 

Different approaches to the NLO calculations of cross sections for hadronic processes 
involving resolved photons exist in literature, see discussion in ||^, |18] and |2^. Here 



we discuss how the NLO QGD calculations, based on the DGLAP approach, are being 
performed for the DIG cross section (fig. [l|,left), 

IP iX, (1) 

where the final photon is produced with large transverse momentum, px ^ ^qcd- 

The Born level contribution to the cross section for the DIG process (|1|), i.e. the 
lowest order in the strong coupling term, arises from the Gompton process on the 
quark (fig. |I|,right ): 

IQ IQ- (2) 

It gives the [Oem] order contribution to the partonic cross section. At the same 
order it contributes to the hadronic cross section for the process 7p — > 7X. The Parton 
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Figure 1. The Deep Inelastic Compton process (left) and Born process (right). 




Figure 2. Examples of single resolved 7 processes: the resolved initial photon (left) 
and the resolved final photon (center) . An example of a double resolved photon process 
(right). 

Model (PM) prediction for the DIG process (|lD, which apphes for xt = '^Pt/ VS ~ 0{1), 
rehes solely on the Born contribution @ namely: 

^^7P-7^ = E / dxpq^p{xp)da^''^^\ (3) 

where qp is the quark density in the proton and da"'^^'^^ [da"''^^'^'^) stands for the 
hadronic (partonic) cross section. In the QCD improved PM the cross section is given 
by (§), however with scale dependent quark densities. For semihard processes, where 
<^ 1, the prediction based on the process (||) only is not a sufficient approximation, 
and one should also consider the contributions corresponding to the coUinear showers, 
involving hadronic-like interactions of the photon(s). There are two classes of such 
contributions: single resolved photon processes with resolved initial or final photon, 
and double resolved photon processes with both the initial and the final photon resolved 
(figs. H). They correspond to partonic cross sections of orders [aemtts] (single resolved) 
and [ctg] (double resolved). If one takes into account that partonic densities in the 
photon and the parton fragmentation into the photon are of order ~ ctem, then the 
contributions to the hadronic cross section from these resolved photon processes are 
al^as and a^^af, respectively. Both single and double resolved photon contributions 
are included in the standard LL QCD analyses of the DIG process 0, |^, ^. To obtain the 
NLO QCD predictions for the process (^ the corrections to the lowest order process 
(0) have to be calculated leading to terms of order a^^as [^, |, ^ (fig. In 
these c^em'^s contributions there are collinear singularities to be subtracted and shifted 
into corresponding quark densities or fragmentation functions. This way the single 
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resolved photon contribution appears in the calculation of the corrections to the Born 
process. It is worth noticing that in the NLO expression for the cross section there are 
no coUinear singularities which would lead to the double resolved photon contributions. 
It indicates that taking into account [af] subprocesses, associated with both the initial 
and final photons resolved, goes beyond the accuracy of the NLO calculation. This will 
be consistent within the NNLO approach, where correction to the Born term and 
correction to the single resolved terms should be included, all giving the same o^em'^s 
order contribution to the hadronic cross sections. 

The other set of diagrams is considered by some authors IH-H and H], H 



in the NLO approach to DIG process (|T]). This approach, which we will call "l/a^" 
approach, is motivated by large logarithms of in the existing already in the 
PM. By expressing In^Q^ / Aq^^j) as ~ l/cts one treats the parton densities in photon 
(and parton fragmentation into the photon) as proportional to aem/c^s (see e.g. 0- 
|],|15[-|1^ and [0, ^ ). By applying this method to the DIG process, we see 



that the single resolved photon contribution to the cross section for the process jp — > 
becomes of the same order as the Born term. The same is also observed for the double 
resolved photon contribution. Namely, we have for the Born, single and double resolved 
photon contributions: 

1 ® FemJ ® 1 = "em^ ® [aemas] ® 1 = a^^, (g) [a^\ (g) = a^^. 

In such counting, the same order contributions to the hadronic cross section are 
given by the direct Born process, single and double resolved photon processes although 
they correspond to quite different final states (observe a lack of the remnant of the 
photon in the direct process). Moreover, they constitute the lowest order (in the strong 
coupling constant) term in the perturbative expansion, actually the zeroth order, so the 
direct dependence of the cross section on the strong coupling constant is absent. Some 
of these terms correspond to the hard processes involving gluons, still there are no terms 
proportional to coupling! 

In the "l/a^" approach the correction to the Born cross section, the single and 
to the double resolved photon contributions are included in the NLO calculation, since 
all of them give terms of the same order, Og^a^, see [|15|-|[1^, [0, 



To summarize, the first approach starts with one basic, direct subprocess as in the 
PM (eq. while the second approach with three different types of subprocesses (as 
in the standard LL calculation). Obviously, some of NNLO terms in the first method 
belong to the NLO terms in the second one. 

In this paper we apply the first type of NLO approach to the DIG process, in 
particular we take into account the following subprocesses: 

• the Born contribution (|^) (fig. |I], left); 

• the finite corrections to the Born diagram (so called K-term) from virtual gluon 
exchange, real gluon emission (fig. ^, left and center), and the process jg qq'j; 

• two types of single resolved photon contributions, with resolved initial and final 
photons (fig. H, left and center). 
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Figure 3. Examples of the virtual gluon (left) and real gluon (center) as corrections 
to the Born contribution. Also a box diagram is shown (right). 



Besides the above full NLO set, we will include two terms of order alrn^l (formally 
from the NNLO set): the double resolved contributions (fig. |^,right) and the direct 
diagram (box) 7(7 — > 7(7 (fig. ^,right), since they were found to be large 

The cross section for the 7p 7X scattering has the following form: 



The first term is the K-term, and the second one stands for the sum over all other 
contributions including the Born term. The /q/^ U'b/p) is a a (6)-parton distribution in 
the photon (proton) while the -D^/c is a c-parton fragmentation function. For the direct 
initial (final) photon, where = 7 (c = 7), we take /^/^ = 5{x^ — 1) (-D^/c = 5{z — 1)) 
(the Born contribution is obtained for a = 'j, b = q and = 7). The variables x^, x and 
z stand for the fraction of the initial photon, proton, and c-parton momenta taken by 
the a-parton, 6-parton, and the final photon, respectively. 

3. The isolation 

In order to observe photons originating from a hard subprocess one should reduce 
backgrounds, mainly from vr^'s and 7's radiated from final state hadrons. To achieve 
this, isolation cuts on the observed photon are introduced in experimental analyses. 
The isolation cuts are defined by demanding that the sum of transverse hadronic energy 
within a cone of radius R around the final photon, where the radius R is defined in the 
rapidity and azimuthal angle space, should be smaller than the final photon transverse 
energy multiplied by a small parameter e: J^hadrons ^Th < ^Et-j 

The simplest way to calculate the differential cross section for an isolated photon, 
daisoi, is to calculate the difference of a non-isolated differential cross section, danon-isoi, 
and a subtraction term, which corresponds to cuts opposite to the isolation cuts 
da sub [EH-O ddisoi = danon-isoi - da sub- 

II Some aspects of the QCD calculation of the isolated photon production are discussed in fsil, |3^. 




abc 



•(4) 
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Note that in practice the isolation cuts are imposed only when calculating the K- 
term, and the contributions involving fragmentation function (resolved final photon). In 
calculation of the subtraction term for the i^-term we applied a small-5 approximation, 
29| , |T5[ 15} Presently we are working on calculations of the cross sections for the 7 



see 



and 7+jet photoproduction using for a comparison the space slicing method |£0[, as in 

r7i,En,Ei. 



4. The results and comparison with data 



We consider the production of photons with large transverse momentum, pt S> ^qcd, 
in the ep scattering, ep — > e'jX, at the HERA collider using the equivalent photon 
( Williams- Weizsacker) approximation |]33| : 



G,/,{y)da^P~'^''dy, 



(5) 



where y is (in the laboratory frame) a fraction of the initial electron energy taken by 
the photon 13^ : 

I + (1 - yf 01^^(1 -y)^ 



-{- 



271 ' y 

with irte being the electron mass. We assume 



m'ly'^ 



y 



2 2 

m^y 



Q 



2 

max 



)}, (6) 



equal to 1 GeV^, as in the recent 
photoproduction measurements at the HERA collider. We neglect the large pt photon 



emission from the electron 35 



The results for the non-isolated and isolated photon cross sections are obtained in 
NLO accuracy with additional NNLO terms, as discussed above. We take the HERA 
collider energies: Ee=27.5 GeV and Ep=820 GeV [|l^, and the pt range of the final 
photon between 5 and 20 GeV {xt from 0.03 to 0.13). The MS scheme with a hard 
(renormalization, factorization) scale Q equal pt (also Q = pt/2 and 2pt) was applied. 
We assume the number of active (massless) flavors to be Nf=A (and for comparison also 
Nf=3 and 5). The two-loop coupling constant is used with Aqc'£)=0.365, 0.320 and 
0.220 GeV for Nf=3, 4 and 5, respectively, as fitted by us to the experimental value of 
as{Mz) = 0.1177 H. 



We use the GRV parametrizations of the proton structure function (NLO and 
LO) the photon structure function (NLO and LO) [^, and the fragmentation 
function (NLO) For comparison other parametrizations were also used: DO @], 

ACFGP i, CTEQ |0|, MRST [|| and GS [H. As the reference we take the GRV 
NLO set of parton distributions |^-||39l, Nf = 4, Aqcd = 320 GeV and Q = px- 



4.1. Non-isolated versus isolated photon cross section 

We have studied the pt distribution for the produced final photon without any cut 
and found that it decreases by three orders of magnitude when p^ increases from 4 

% This small S approximation seems to be an accurate analytic technique for including isolation effects 
in NLO calculations (also for i? = 1), see discussion in Q. 
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GeV to 20 GeV (not shown). Obviously the most important contribution is coming 
from the lowest px region, where the resolved photon processes dominate. The total 
NLO cross section integrated over pt range from 5 to 10 GeV, is equal to 226 pb, with 
individual contributions equal to: Born = 36.3%, single resolved = 35.1%, double 
resolved = 18.7%, box = 6.2%, K-term=3.9%, so the single resolved photon processes 
give a contribution comparable to the Born term. Also the double resolved photon 
processes are important. The direct box diagram {jg 7(7) gives 17% of the Born 
ill ~^ 11) contribution. Such relatively large contribution is partially due to large 
gluonic content of the proton at small Xp. 

The differential cross section for the final photon rapidity, da/dr/^, for the non- 
isolated photon and for the photon isolated with various cones (various e, R) was 
studied. The isolation cut suppresses the cross section significantly in the whole rapidity 
range(not shown). For e=0.1 and R = 1 the suppression is 17-23% at rapidities 
-1.5 < r^^ < 4 0. 

As expected, the cross section for fragmentation processes is strongly suppressed: 
after isolation it is lowered by a factor of 5. At the same time the QCD corrections to 
the Born diagram increase significantly, i.e. the contribution to the subtraction cross 
section, dagub, due to this corrections is negative. The subtraction cross section, being a 
sum of negative QCD corrections and fragmentation contributions, is of course positive. 



4-2. Other experimental cuts 

In order to compare the results with data we fix R=l and e=0.1, which are the standard 
values used in both theoretical and experimental analyses, and consider other cuts 
imposed by the ZEUS group |jl2[. Two types of the final state were measured in the 



ZEUS experiment: 1) an isolated photon with —0.7 < rj'^ < 0.9 and 5 < < 10 GeV; 
2) an isolated photon plus jet with the photon rapidity and transverse momentum as 
above, the jet rapidity in the range —1.5 < 77-''^* < 1.8, and the jet transverse momentum 
Pt* > 5 GeV. Here we compare our NLO predictions with the ZEUS data from the first 
type of measurements [|T2|. More results can be found in and . 



We have found (see [|T9[) that the cross section for a production of final 7 is strongly 
reduced, by 30-85%, in the positive rapidity region due to the limited energy range, 
0.2 < y < 0.9 (see also |^T], ^ for a similar conclusion). At negative rapidities the 



change due to the y-cut is weaker: 5-10% at —1.2 < rj^ < —0.4 and 10-30% at other 
negative rapidities. Note however, that the Born term is reduced 3.5 times. 

The role of various experimental cuts is illustrated in fig. ^ for the x^ distribution. 
The small x^ contributions are strongly, by two orders of magnitude, diminished by 
the photon rapidity cut. This shows that measurements at the central rj'^ region 
(—0.7 < 7]"' < 0.9) are not too sensitive to the small values in the photon, and 
can not be used presently for constraining eg. the gluon density in the photon. 

+ The positive rapidity is pointed in the proton direction. 
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Figure 4. The cross section in bins of the length 0.1. The results for non- isolated 
7 are shown (dashed line). The solid line represents results for isolated 7 with e = 0.1 
and R — 1. Results with additional cuts in the isolated 7 cross section are shown with: 
dotted line (0.2 <y< 0.9) and dot-dashed line (0.2 <y< 0.9, -0.7 < 7/.^ < 0.9). 



ep — > e^fX ep — > e^^fX 




Figure 5. The results for isolated 7 for: TV/ = 3 (dashed lines), 4 (solid lines) and 
5 (dotted hnes), compared to the ZEUS data The da/dpr as a function of the 
photon transverse momentum (left) and du/drj^ as a function of the photon rapidity 
(right); the result without the box contribution is also shown for Nf = 4 (dot-dashed). 



4-.3. The comparison with data 

In fig. ^ (left) tlie comparison is made witli ZEUS data for the pt distribution for isolated 
7 for various Nf. A satisfactory agreement is obtained for Nf = 4 (and 5). Note large 
difference between the results for Nf=4 and 3 due to the fourth power of electric charge. 
The rapidity distribution is shown in fig. |^(right), where a good description of the data 
is obtained for iVj=4 and 5 in the rapidity region 0.1 < rj^ < 0.9. For —0.7 < rj^ < 0.1 
our predictions lie mostly below the experimental points. This disagreement between 
predicted and measured cross sections is observed also for other theoretical calculations 
(LG,FGH) and for Monte Carlo simulations [|12|. In fig. ^(right) we present separately 
an effect due to the box subprocess (for Nf = 4). It is clear that the box term enhances 




Figure 6. The da/drj-y for isolated 7 compared to tlie ZEUS data [12|. Three different 
NLO parton distributions for 7: ACFGP [|| (dotted hne), GRV 13?] (sohd Hne) and 
GS II (dashed hne). Q = pt (left) and Q = 2pT (right). 



considerably the cross section in the measured rapidity region (by ~ 10%). The double 
resolved photon contribution is also sizeable, although roughly two times smaller than 
the box one. Both these [a^] contributions improve description of the data. 

The predictions obtained using three different NLO parton densities in the photon 
(ACFGP §, GRV [0 and GS []^) are presented for Nj = 4 in fig. | for Q = pr 
(left) and for Q = 2pT (right) together with the ZEUS data |12|. The results based 
on ACFGP and GRV parametrizations differ by less than 4% at rapidities rj^ < 1 (at 
higher rj^ the difference is bigger), and both give good description of the data in the 
rapidity range 0.1 < r/^ < 0.9 (for Q = pr and Q = 2pt)- For —0.7 < r]^ < 0.1 none of 
the predictions is in agreement with the measured cross section. For Q=Pt (fig- Heft) 
the GS distribution leads to results considerably below ones obtained using ACFGP 
and GRV densities, especially in the rapidity region from roughly -1 to 1. This is due 
to a different treatment of the charm quark in the photon, namely in the GS approach 
the charm quark is absent for below 50 GeV^ - contrary to GRV and ACFGP 
parametrizations where the charm threshold occurs at lower Q^. All the considered 
parton distributions give similar description of the data when the scale is changed to 
Q = 2pT, since then Q"^ is always above 50 GeV^. The (data-theory) /theory for the 
same cross section for the reference set (GRV) of parton parametrizations is presented in 
fig. 1^. In fig. H our predictions are compared to the ZEUS data divided into three ranges 
of y. Clearly the discussed above discrepancy between the data and the predictions for 
rj"' < 0.1 is coming mainly from the low y region, 0.2 < y < 0.32. In the high y region, 
0.5 < y < 0.9, a good agreement is obtained for a whole measured rapidity region. This 
is not the case of LG and FGH results, which are in disagreement even for a large y 
(mainly for rapidities above 0.1, see below). 

We have also studied the dependence of our results on the choice of the parton 
distributions in the proton and parton fragmentation into the photon (not shown), and 
a small sensitivity was found. Only at minimal [t]^ < ~1) and maximal (4 < rj^) rapidity 
values this difference is larger, being at a level of 3.5 — 8%. 
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Figure 7. (Data-theory) /theory is plotted for GRV parton parametrizations for the 
initial and final photons, and for proton, the ZEUS data from [p2|. 
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Figure 8. The results for three ranges of y: 0.2 < y < 0.32, 0.32 < y < 0.5 and 



0.5 < y < 0.9, compared to the ZEUS data |12 



5. The theoretical uncertainties and comparison with other NLO analyses 

In order to estimate the contribution due to missing higher order terms, we have studied 
the influence of the choice of the Q scale for the rj^ distribution. Some of results can 
be found in flg. ^ for the GRV and ACFGP parton parametrizations. Around the 
maximum of the cross section at rapidities — 1 < < changing Q scale from pt/2 
to 2pT leads to differences 4-6%. This small sensitivity of the results to the change of 
the scale is important since it indicates that the contribution from neglected NNLO and 
higher order terms is not signiflcant. Note that individual contributions are strongly 
dependent on the choice of Q, e.g. results for the single resolved processes vary by ±10- 
20% at rapidities rj^ < I. Results are much more stable only when the sum of resolved 
processes and QCD corrections is considered. These results leads to expectation that 
our prediction should not differ too much from results based on larger set of diagrams. 

As we discussed in Sec. ^ our NLO calculation of the DIG process differs from the 
"l/as"-type NLO analysis presented in ref. [|l5]-[|l3 and I^T], |2^, by set of diagrams 
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Figure 9. The comparison of tiie data for the px distribution (from |T^) with 
PYTHIA and HERWIG Monte Carlo predictions and with our results and LG ones 
based on GRV parametrizations (left). The same data in comparison with the FGH 
results (left, from ||2l|). 



included in the calculation. We do not take into account corrections to the single 
and double resolved processes, which are beyond the NLO accuracy in our approach. 
On the other hand, we include the box diagram neglected in |[T^-[0, which however 
is taken into account in FGH analysis (The double resolved subprocesses are 

included in all mentioned analyses.) 

We compare now our results and the results of the LG calculation (using A^/=4 
and Q = px) for the isolated final photon (-R=l, e=0.1) in the kinematical range as in 
the ZEUS analysis [|12[ (i.e. for -0.7 < r/T < 0.9 and 0.2 < ?/ < 0.9), see fig. |(left). In 
fig. ^(right) a comparison is made for the FGH results ||2T[ and the same data. The LG 
predictions (and FGH results) for da/dpr (with GRV parton parametrization) are about 
20% higher than ours in the presented range of transverse momentum, 4 < p-r < 20 
GeV. For da/drp cross section (with 5 < < 10 GeV) the biggest differences between 
our and LG predictions are at large y range, what can be seen in fig. |10| (left). For y 
range limited to low values only, 0.2 < y < 0.32, the LG cross section is higher than ours 
by up to 20% at positive 77^, while at negative rj^ it is lower by up to 10%. For large y 
values, 0.5 < y < 0.9, where our predictions agree with data, the LG results are higher 
than ours by up to 80% (at rj^ = 0.9). Not only the LG predictions are too high as 
compared to the data in the forward direction, similar effect especially for large y is seen 
in the fig. lO(right), where the FGH results are compared with data. Our predictions 
as we discussed above are close to the data for i]^ above 0.1, see figs. ^ | and |10|(left). 



6. Summary 

Results of the NLO calculation, with additional NNLO contributions from double 
resolved photon processes and box diagram, for the isolated 7 production in the DIG 
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subdivisions of photon energies 




Figure 10. The comparison with data for the rapidity distribution of the final photon 
three ranges of y (as in fig. 7 ) of PHYTIA prediction, our results and LG results (for 
GRV), (left, (from jlj])) and similar comparison between data and FGH results, from 
P (right). 



process at HERA are presented The role of the kinematical cuts used in the ZEUS 
measurement are studied in detail. 

The results obtained using GRV parametrizations agree with the data in shape and 
normalization for px distribution. For t]"' distribution a good description of the data 
is obtained for 77'^ > 0.1, while for 77"^ < 0.1 the data usually lie above the predictions. 
This discrepancy arises mainly from the low y region, 0.2 < y < 0.32. The beyond NLO 
terms, especially a box contribution, improve the description of the data. 

We have studied the theoretical uncertainty of results due to the choice of the 
renormalization/factorization scale: Q = pt/2,pt,2pt- At high rapidities rj^ > 3, 
where the cross section is small, this uncertainty is 10-30%. In a wide range of rapidities, 
—2<ri^< 2, the dependence on the Q scale is small, below 6%. Since we include some 
NNLO diagrams in our NLO calculation, this stability of the predictions versus the 
change of the scale is especially important. The week dependence on the Q scale, and 
not large differences between LL and NLO predictions (below 20%) allows to conclude 
that theoretical uncertainties of our NLO calculations for an isolated photon production 
in the DIG process at HERA are relatively small. 

We compared our results with the "l/ag" NLO calculations by LG and FGH, which 
are based on different set of subprocesses. The cross section da/dpT obtained by LG is 
about 20% higher than ours (for GRV photonic parton distributions), FGH prediction 
is closer to ours than the LG one, it lies between KZ and LG curves. For the cross 
section da/drj^ the difference between our results and LG/FGH ones is up to 35% at 
rj^ = 0.9. The highest differences are present for high y values only, 0.5 < y < 0.9, 

* Our fortran code is available upon request from azem@fuw.edu.pl. 
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where on the other hand our predictions are in agreement with the data. At low y 
range, 0.2 < y < 0.32, differences between our calculation and calculations done by LG 
and FGH are smaller and none of them describe the data well for rapidities below 0.1. 
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